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TECHNICAL PUBLICATION 


PERFORMANCE THEORY OF DIAGONAL CONDUCTING WALL 
MAGNETOHYDRODYNAMIC ACCELERATORS 

1. INTRODUCTION 


Historically, interest in magnetohydrodynamic (MHD) devices has centered on their use as electri- 
cal generators in commercial central power plants and mobile burst power systems. The primary attraction 
for central power is associated with the attainment of higher peak cycle temperatures, which point to 
significant improvements in overall plant efficiency. The attraction for mobile pulse power, particularly 
airborne military needs, derives from their intrinsic high power-density characteristics. It should be recog- 
nized that equally important though less noted attractions exist for accelerator configurations, as well. 

Two major identifiable uses for crossed-field MHD accelerators are (1) propulsive devices and 
(2) hypersonic aerodynamic test facilities. Litchford et al. discuss the operational attributes that are of 
particular significance to these applications and provide an indepth historical perspective of their techno- 
logical development. 1 

From a fundamental point of view, the essential argument favoring the utilization of a Lorentz force 
acceleration mechanism is the ability to avoid the inherent physical limitations encountered with pure 
thermal approaches (viz, material thermal limits and ionization/dissociation losses). Simply put, it is more 
effective to transfer electrical energy into directed kinetic energy rather than first degrading it to thermal 
energy. Crossed-field acceleration is of special interest in these cases because of its unique capacity for 
processing large amounts of power under conditions of high mass throughput. 

Past emphasis on generator configurations has led to a substantial theoretical and experimental 
basis for understanding their performance and operational nuances. However, a similar level of under- 
standing has yet to be acquired for accelerator configurations despite a strong growing interest in their 
potential aerospace applications. 

The standard theoretical approach for describing the interaction of an electrically conducting gas 
with applied electric and magnetic fields relies on the application of a Cowling-Schlueter-type generalized 
Ohm’s law for a partially ionized, electrically neutral gas. Powers et al., for instance, used a cross-plane 
averaged form of the generalized Ohm’s law and developed a graphically based methodology for describ- 
ing the general performance characteristics of MHD accelerators and generators. 2 

This classical model has proved to be extremely useful in illuminating the basic operational behav- 
ior of MHD devices, but its practical utility is hampered by a critical limitation. That is, the theory assumes 


a constant near-electrode voltage drop, whereas experience has shown that the magnitude of the near- 
electrode voltage drop exhibits a significant load-current dependency. 

Various theoretical refinements of this nature have been forthcoming but solely within the context 
of generator configurations. Wu, for example, introduced an effective voltage drop parameter, accounted 
for its load current dependency, and examined the resulting impact on diagonal conducting wall (DCW) 
generator performance . 3 The MHD research group at the Max-Planck-Institut fuer Plasmaphysik made a 
similar attack on the problem using a cross-plane averaged Ohm’s law, specialized for a DCW generator 
configuration . 4 

It is widely recognized, of course, that the current and electric field structure in MHD devices is 
inherently three-dimensional and that spatial separation of physical processes is not warranted in general. 
As a matter of recourse, one could mount a brute force attack on the problem and perform a complete three- 
dimensional numerical analysis of the flow and electrical structure in the duct, assuming adequate comput- 
ing resources are available. Alternatively, one could invoke an infinite segmentation assumption (i.e., the 
streamwise variation in electrical properties is small in comparison to transverse variations) and couple an 
approximate cross-plane electrical model with a three-dimensional flow analysis as exemplified for 
generator configurations by Bityurin et al. and Ahluwalia et al . 5 ’ 6 Each approach has merits and draw- 
backs, and experience has shown that each level of approximation is useful in fulfilling certain design and 
analysis needs. 

The purpose of this Technical Publication is to reexamine the classical theory and extend it for 
cross-field accelerator configurations with inclusion of a current-dependent near-electrode voltage drop 
model. This particular refinement of the classical theory, while generally recognized as a straightforward 
theoretical extension, has not explicitly appeared in the literature, and it is believed that the resulting 
development can yield practical insight into the basic operational characteristics of these important 
devices. Moreover, it is hoped that this work may serve as a convenient and compact resource for future 
design practitioners. 
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2. CROSSED-FIELD ACCELERATOR CONFIGURATIONS 


Various linear crossed-field MHD accelerator configurations are conceivable (fig. 1) although the 
optimal configuration depends on the ultimate application need. The Hall configuration, for example, is 
generally more effective for low-density flows (i.e., very high Hall parameter condition) whereas the 
segmented Faraday configuration yields superior performance for high-density flows. A significant disad- 
vantage of the Faraday configuration, however, is the separate power conditioning that is required for each 
electrode pair, which leads to a heavy, complex, and expensive system. 

Alternative two-terminal loading schemes have been proposed to avoid the complications of 
multiterminal connections while also attempting to compensate for the effects associated with a finite Hall 
parameter. De Montardy, for instance, suggested the series connected scheme in which a segmented 
Faraday channel is externally diagonalized. 7 Dicks later extended this approach to DCW configuration in 
which slanted window frame-like electrode elements are stacked with thin insulators to form a complete 
channel. 8 ’ 9 The DCW configuration not only simplifies fabrication and improves strength but provides 
superior performance to the series connected device by allowing current to flow to the sidewalls. 3 

Hence, theoretical development is focused on the DCW accelerator configuration, since it has the 
greatest practical relevance. The performance characteristics of the alternative crossed-field configura- 
tions can be subsequently deduced as special or extreme cases of the DCW configuration. 
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3. GENERALIZED OHM’S LAW 


The local current density, j, in a partially ionized, electrically neutral gas may be determined from 
a Cowling-Schlueter-type generalized Ohm’s law. 10-14 For application to typical MHD devices, it is useful 
to adopt the simplified form derived by Rosa, which includes ion slip but neglects electron pressure 
gradients: 15 


j = aE'-^jxB + ^r(jxB)xB . (!) 

B B 2 

Here, o is the scalar electrical conductivity; E' is the generalized electric field intensity relative to axes 
moving with gas velocity, u; B is the magnetic field strength; and fl e and /3 ; are the electron and ion Hall 
parameters, respectively. The Hall effect and ion slip are introduced by the second and third terms on the 
right-hand side of equation (1). 

As a practical matter, we define the generalized electric field, E', in terms of the electric field in the 
laboratory frame, E; the motional emf, uxB; and an equivalent electric field, E d , which accounts for the 
near-electrode voltage drop, V d (i.e., E d =V d /h), inside the boundary layer, 

E' = E + uxB-E^ . ( 2 ) 

Note that the equivalent electric field is applied only in a direction transverse to the magnetic field vector 
and the streamwise velocity vector. It is always directed in opposition to the transverse current, j . 

For the special case in which no current flows along the magnetic field lines, ( /_=()), it is readily 
shown that equation (1) reduces to the compact form 


Q 

j = £ E' jxB , 

j B 

where we have introduced two new parameters: 


(3) 


i + A?A 


a 

1 + b 


and 


Pe _ Pe 

1 + A'A 1 + b 


(4) 


( 5 ) 
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In many cases, the ion slip factor, b=f5 e /3 i , is negligible, and further simplification is obtained as Zetland 
Cl->P e . 


3.1 Mean Ohm’s Law 


Many of the interesting phenomena that occur in high-interaction MHD devices, such as velocity 
overshoots in the sidewall boundary layer, flow asymmetries, and generation of intense secondary flows in 
the cross plane, are three-dimensional in nature. Therefore, high-fidelity performance predictions ulti- 
mately require three-dimensional numerical analyses of the combined flow and electrical structure. On the 
other hand, simplified approaches based on spatial decoupling or reduced spatial dimensions can also yield 
useful results if all critical physical effects are accounted for in an appropriate fashion. 

Exact theoretical treatment requires consideration of the fundamental equations of MHD 
and Maxwell’s electromagnetic equations, thereby leading to a system of coupled differential equations. 
Here, we invoke the infinite segmentation assumption and assume that all parameters in the generalized 
Ohm’s law can be effectively represented by their cross-plane averaged values and are dependent on the 
x-coordinate, only. 

The dimensionless voltage drop parameter, A =E d /uB=V d /uBh, which incorporates all voltage losses 
associated with the electrode boundary layer, may be introduced into the generalized Ohm’s law to obtain 
the component relations, 


Equations (6) and (7) may also be combined to obtain the inverse expressions for j x and j in terms of the 
local electric field components, 


jx~^E x Qjy 


(6) 


and 


jy = 'LE y -'LuB(l + A) + Qj x . 


(7) 



E Y — E v + 


y YjUB{\ + A) 

1 + Q 2 


( 8 ) 


and 



y2kB(l + A) . 

1 + Q 2 


(9) 
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4. PERFORMANCE MODEL 


4.1 Diagonal Loading Constraint 

Completion of the electrical model requires the consideration of external loading conditions. In 
the case of a diagonally linked device (fig. 2), the electric field is forced to align perpendicularly to the 
diagonal short so that 


Ey jE x = tan0= (p , (10) 

where 6>is the electric field angle and (p is the electric field direction. It follows that the diagonalization or 
wall angle is given by 6 w =n/2-6. The sign of depends on the mode of operation. For an accelerator, uxB 
opposes j , 0<0 <k/2, and (p> 0. For a generator, uxB aligns with j , -7t/2<0<(), and (p<0. 

The two-terminal load current, /, for an MHD device with diagonally linked electrode pairs is given by 


l = \ Af > •rfA / = j "A / , 

where the integration is over the entire slanted area, A ,-. In component form, 

I = ( j x + tan °jy) A = (jx + Wy ) A • 


( 11 ) 


( 12 ) 


4.2 Electrical Parameter Relationships 

A complete set of equations now exists for determining the cross-plane electrical characteristics. In 
general, either I or E x may be specified, allowing the remaining unknown electrical parameters to be 
deduced from equations (6), (7), (10), and (12) in conjunction with appropriate material functions. 

The resulting performance relationships for DCW accelerators are summarized in table 1 using E x 
and I as independent variables. The governing relations for generator performance are recovered when 
both (p and A are less than zero and I=-E X L/R, where R is the load resistance and L is the channel length 
over which the load is applied. Note also that the governing relations for a linear Hall channel configura- 
tion are recovered in the extreme case (p=0. 
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4.3 Mean Power Parameters 


Ultimately, we seek to deduce suitably averaged power parameters from the resulting cross-plane 
electrical model. The quantities j-E and (jxB) r , for instance, may be used to determine the power density, 
push power, dissipated power, and electrical efficiency of a device. 

4.3.1 Power Density 

The electrical power density at any cross section of the duct is defined by 

P=i- E = ^Ji!+d • (13) 


Eliminating E using equation (10) yields 


P = 


Aj Af 


but given A f 



, the working form becomes 


P = 


IEx 

A 


(14) 


(15) 


4.3.2 Push Power 

The streamwise Lorentz body force component at any cross section is defined as 

(jxB ) x = jy B , (16) 

and the push power associated with this Lorentz body force is given by 

P p = u- (jxB) x = uj y B . (17) 

4.3.3 Dissipated Power 

The power density of an MHD device must exceed the Lorentz force push power since the internal 
resistance to current flow leads to ohmic heating. This can be demonstrated by taking the dot product of j 
with the generalized Ohm’s law, as defined by equation (3), which yields 
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( 18 ) 


j j = Z[j-E + j-(uxB)-j Ej-^j- (jxB) . 

D 

Because /'_=(), the last term is identically zero, and we obtain 

-jjr = j'E-u- jxB- yE d =P-P p -j y uBA , ( 19 ) 

where jj/X is the power dissipated in the core flow region, and j -E d =j E d =j uBA represents the power 
dissipated in the electrical boundary layers. Therefore, the power dissipated over the entire cross plane is 
given by 


Pd 


jj 

X 


+ j y uBA=P-Pp 


( 20 ) 


4.3.4 Electrical Efficiency 

The electrical efficiency of an accelerator is simply the ratio of the push power to the applied 

power: 


P p u ■ (j x B) t _ ujyB 
P ' j-E IE X /A 


( 21 ) 


It follows that the generator efficiency is defined as the reciprocal of the accelerator efficiency, r] ? =l/rj fl . 

4.3.5 Power Relationships 

By combining the results in table 1 with the relationships above, it is possible to express the mean 
power parameters in terms of the independent variables E x or 7. These results are summarized in table 2. 


4.4 Generalized Vector Diagram 

It is generally instructive to examine the qualitative behavior of the internal electric fields and 
current densities under the assumption that MHD effects have only a slight effect on flow properties. In this 
case, the generalized Ohm’s law can be used to construct representative vector diagrams for diagonalized 
devices, as shown in figure 3 for fixed values of (p. 

To construct these diagrams, the generalized Ohm’s law defined by equation (3) is put into the form 



i*» 

X B 


= E' 


( 22 ) 
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where we immediately recognize that the two terms on the left-hand side of equation (22) must be orthogo- 
nal. Thus, each of these terms represents one side of a right triangle with the hypotenuse, E', inclined at an 
angle of tan -1 Q. The generalized electric field must simultaneously satisfy the vector relation defined by 
equation (2). 

Figures 3(a)-(c) illustrate the load current dependence of a diagonalized generator ((p< 0). In the 
open-circuit case (a), there is no current flow to the load, and the internal current density vector must align 
with the diagonal link, which defines planes of constant electric potential. Thus, the current is closed 
through the diagonal short, and E exactly compensates for the combined uxB and E rf induced potentials. 
Note that j is directed upstream under these conditions. As the load current is increased, the current 
density vector rotates out of the plane of the diagonal link until it aligns with the uxB vector. At this point, 
j x = 0 and we arrive at the Hall current neutralized condition (b). If the load current is increased further, the 
current vector continues its rotation, and the direction of Hall current flow is reversed. Ultimately, the 
current density vector becomes perpendicular to the plane of the diagonal link as the electric field goes to 
zero at the short-circuit condition (c). 

Figures 3(d)— (f) illustrate the load current dependence of a diagonalized accelerator (q» 0). When 
the applied axial electric field is zero and there is no load current, the open-circuit accelerator loading 
condition is obtained (d). The internal current density vector corresponds to the generator short-circuit 
case, under these unique circumstances. Note that j < 0 (braking regime) and that the current flow is closed 
through the diagonal short. With increasing applied axial electric field, the current density vector rotates 
out of the plane of the diagonal link. The device remains in a braking regime, however, until the applied 
field becomes large enough to force j to zero at which point we attain the zero-efficiency accelerator 
condition (e). Further increases in applied current lead to continued rotation of the current density vector 
and increasing values of j . Eventually, one reaches a j x neutralized- accelerator condition (f). Increasing 
the applied current beyond this point leads to a reversal in the direction of axial current flow while the 
magnitude of the current density vector continues to grow. 

4.5 Current-Dependent Voltage Drop 

The previously developed performance relationships hold for any fixed value of A. However, 
extensive experimental research with MHD generators has clearly established that the near-electrode volt- 
age drop is approximately linearly dependent on the transverse current density. 16-18 That is, any increase 
in load current leads to a proportionate increase in voltage drop. This behavior is assumed to extend to 
accelerator operation, as well. 

The linear relationship between V d and j permits the construction of a simple two-parameter model, 


V* 


Vd ~Vd 

I" -I' 


{i-n+Vd 


(23) 


where the prime and double-prime superscripts indicate two distinct loading conditions. It is generally 
more convenient to express this relation directly in terms of A: 
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(24) 


A = 


A -A 


(/-/') + A . 


Selection of appropriate reference load conditions depends on the mode of operation. For a genera- 
tor, it is most natural to adopt the open-circuit (A'=A 0 ) and short-circuit (A"=A V ) conditions on the extreme 
ends of the operational load line. For an accelerator, it is convenient to adopt the zero-power condition, 
which corresponds to a short-circuit generator (A'=A_=A^), and the j neutralized condition (A"=A /; =0). 
The selected reference conditions are summarized in table 3. The implications of this model are now 
examined for the generator and accelerator modes. 


4.5.1 Generator Mode 


At the short-circuit generator condition, E x = 0, 1"=I S , and A"=A v . From table 1, we therefore deduce 



'LuBA(l + A s )(Q-(p) 
1 + Q 2 


(25) 


Elimination of I" using equation (24) immediately yields an expression for A in terms of / for a generator: 


A = A 0 


(i+q 2 )K-a.,) 1 l 

(Q-<p)(l + A ? ) ZuBA 


(26) 


This relationship may be used to eliminate A from the performance relations of tables 1 and 2. 

4.5.2 Accelerator Mode 

The zero-power open-circuit accelerator condition corresponds to the short-circuit generator con- 
dition with cp> 0. In this case, E = 0, /'=/_, and A'=A_=A V , and we deduce that 


r 


= i 


z 


ZuBA(l + A s )(n-(p) 
1 + Q 2 


(27) 


The point of transverse current-density neutralization corresponds to a zero-efficient accelerator where 
j y = 0, 1"=I n , and A"=A w =0. This implies 



YjiiBA 
Cl + (p 


(28) 


10 


Eliminating /' and I" in equation (24) gives an expression for A in terms of / for an accelerator, 


(Q + (p)(l + Q 2 

) ' -i 

'hi BA ' 

(q 

2 -„ 2 ) 

|(1 + A,) 

(l + Q 2 ) 

-(n 2 V) 

|(1 + A,) 


This relationship may be used to eliminate A in the performance relations of tables 1 and 2. 


5. PERFORMANCE DIAGRAMS 


It is possible to construct an informative performance map for DCW MHD devices following the 
general methodology first outlined by Powers et al., who assumed a constant voltage drop in the electrode 
sheath layer. 2 The significant modifications associated with a variable effective voltage drop model were 
later noted and described by Wu for the DCW generator case. 3 Here, Wu’s analysis is explicitly extended 
for DCW accelerator operation, as well. 

5.1 Dimensionless Performance Parameters 

As a first step, the governing performance relationships are simplified by defining the following 
dimensionless quantities: 


-j=-L , e=*.1=-L- .p—^2 ■ 

EuB uB 2 uBA Eu~B 2 


(30) 


Substitution into the component form of the generalized Ohm’s law then yields 


Jx 1 + Q 2 


E x +QE y +£l(l + A) 


(31) 


and 


Jy 1 + Q 2 


£IE X +E y -(1 + A) 


(32) 


It also follows that the electric field components can be stated explicitly in terms of current density: 


E x = j x +Qj y 


(33) 


and 


Ey - ~£lj x + jy+(l + A) . 


(34) 


An expression for the dimensionless load current is obtained from the substitution of equation (30) into 
equation (12): 


I = j x +(pj y , 


(35) 
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and the dimensionless laboratory power density follows directly from the substitution of equation (30) into 
equation (13): 


P = j E = j x E x + j y Ey . 

Eliminating E x and E using equations (33) and (34) yields the more convenient form, 

P = j 2 x+jy+(l + A)j y 


(36) 


(37) 


Furthermore, we may define the dimensionless push power as 

- _ Pp _ 11 .! y P _ jy _ — 

P 'LirB 2 Zu 2 B 2 2 uB J y ’ 


(38) 


and the electrical efficiency may be expressed in terms of dimensionless parameters by combining 
equations (37) and (38): 


P P _ j y _ 1 

P jl + 7v + ( 1 + A )Jy r>8 


(39) 


It is readily apparent that the relationships in tables 1 and 2 can all be put in dimensionless form using 
equation (30); the results are obvious and are not tabulated. The current-dependent effective voltage drop 
effect is accounted for by using either equation (26) or (29) in dimensionless form. 

As a prelude to the construction of the graphical performance map, it is instructive to examine 
performance parameter variations as a function of the dimensionless axial electric field, E r Figure 4, for 
example, summarizes calculated parameter variations assuming representative values for electrical/flow 
conditions, material functions, and wall angle. Starting from the open-circuit generator condition and gradu- 
ally reducing the load impedance, we observe the sequential occurrence of all of the loading conditions 
depicted in figure 3. 

For an open-circuit generator, 7=0 and equation (35) yield the constraint j x 0 =-(pj x 0 . We therefore 
infer that j x o <0 since cp< 0 and j o <0 for the assumed magnetic field direction. In this case, the internal 
current is forced parallel to the wall angle and is completely shorted through the diagonal linkage. As the 
load impedance is decreased below a critical level, however, the total current becomes nonzero, and the 
magnitude of the negatively directed axial electric field begins to fall. 

If the load impedance is decreased all the way to zero, we arrive at the short-circuit generator 
condition corresponding to F v = 0. Here, application of equation (33) yields the constraint 
j x =-Q j y s , which implies j xs > 0 since Q>0 and j < 0. In this case, the internal current flow is forced 
perpendicular to the wall angle, and no current flows through the diagonal linkage. 
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Consideration of these open- and short-circuit constraints clearly implies the existence of a Hall 
current neutralized condition (i.e.,y v =0) on the generator load line. It is also clear that maximum power 
extraction and electrical efficiency must also occur at some point between the open- and short-circuit 
conditions. This point may or may not coincide with the Hall current neutralized condition. Note that the 
magnitude of A is largest at open-circuit loading where the magnitude of j is greatest. 

Simply reversing the wall angle (<p>0) at the generator short-circuit condition yields the open- 
circuit accelerator condition. If power is then gradually applied such that a positive axial electric field is 
impressed on the accelerator, we enter a regime where the transverse voltage is insufficient to overcome 
the induced uxB potential, and j remains negatively directed. This is commonly referred to as the braking 
regime, where all applied power goes into joule heating of the working fluid. Eventually, the applied 
voltage becomes high enough to neutralize the transverse current at which point A— >0. 

Increasing the applied electric field beyond this point yields positive push work and flow accelera- 
tion. The transverse current and total current grow in magnitude, more power is delivered to the device, 
and A steadily increases. At the same time, axial current flow gradually decreases until the j neutralization 
condition is reached. Additional increases in applied field and power beyond this point simply reverse the 
direction of j . It is important to note that optimal accelerator electrical efficiency occurs at an applied 
electric field magnitude much less than that required to neutralize axial current flow. 

5.2 Performance Map 


The operational attributes of MHD devices can be exhibited in the most vivid manner on a perfor- 
mance map in the plane of dimensionless current density and electric field. Following the methodology of 
Powers et al., 2 as it was modified by Wu for current-dependent voltage drop, 3 we complete the square in 
equation (37) to obtain the family of circles defined by 



jy+- 


1 + A 


= P+ 



(40) 


These circles are centered at j x =0,j y =(\+A)/2 with a radius of {P+[(l+A)/2] 2 } 0-. 

Transformation to the E x , E y axes is defined by equations (33) and (34), and the origin of the E x , E y 
axes in the j x ,j plane follows directly from equations (31) and (32): 


j.x,E = 0 


Q(l + A) 
1 + Q 2 


7 (l + A) 

y - E =° 1 + Q 2 ' 


(41) 


The locus of points defining this origin can be expressed independently of Q by combining equations (33) 
and (34) under the condition E X =E =0: 


72 _ , 

J x,E=0+ 


Jy,E= 0 + 2 


l + A 


l + A 


(42) 
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Note that the circles defined by equation (40) reduce to the loci defined by equation (42) at zero-power 
density. Thus, the generator boundary coincides with the locus of the origin of the E v E y axes, and all 
power generating circles lie inside this boundary since P < 0 for a generator. Further insight may be gained 
through the examination of extreme generator loading conditions. 

At generator short circuit, the Hall field goes to zero, and equation (33) gives the constraint 

Because Q>0 and j <0, given the assumed magnetic field orientation, the short-circuit operating domain 
resides in the fourth quadrant of the current density plane. 

At generator open-circuit conditions, the total current falls to zero and equation (35) gives the 
constraint 



Because <p < 0 and j y < 0 for the generator mode, the open-circuit operating domain resides in the third 
quadrant of the current density plane. 

Clearly, the power density vanishes at both open and short circuit conditions, but the critical 
discriminating factor to note is that A varies with current over the entire load range (-l<A o <A s <0). There- 
fore, the radius of the zero-power circle depends on the load. The zero-power loci for short- and open- 
circuit loading follows from equation (40) with P= 0: 

-2 . 

J ws" 1 " 

and 

-2 , 

J x, ca- 
using these expressions, it is possible to construct the short- and open-circuit half-circles. The 
short-circuit half-circle is centered at j x =0, j y =(\+A s )/2 with radius ( 1 +A v )/2, and the open-circuit half- 
circle is centered at j x =Q, j =(l+A 0 )/2 with radius (l+A 0 )/2. These loci are illustrated in figure 5 where the 
load line connects in accordance with the constraints defined by equations (43) and (44). Because A 0 <A V , 
the radius of the open-circuit half-circle is smaller than the radius of the short-circuit half-circle, and the 
center points are shifted. 

The origin of the E x , E y axes falls on the short-circuit half-circle at the points defined by equation 
(41). Inspection of equation (34) indicates that the E axis passes through the origin of the j x ,j axes since 
A=0 at this condition. The various points of interest are indicated in figure 5. 




1 + A, 


1 + A, 


(46) 


Jy,s + ' 


1 + A , 


1 + A, 


(45) 
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Combining the DCW electric field constraint defined by equation (10) with equations (30), (33), 
and (34) yields an expression for the operating load line in the form 


It is now possible to eliminate A using either equation (26) for a generator or equation (29) for an accelera- 
tor. In both cases, one ends up with an equation for a line in the j x , j plane since A is linearly related to j x 
and j . The resulting operating lines are illustrated in figure 5. 

In the generator case, the operating line passes through the short- and open-circuit load conditions, 
as determined by equations (43) and (44), respectively. In the accelerator case, the operating line originates 
from the generator short-circuit loading point on the zero-power half-circle. This line also projects through 
a point on the half-circle defined by the value of cp. 

The P circles for accelerator operation are obtained from equation (40) using equation (29) to 
eliminate A. Note that the accelerator operates in the dissipative braking regime between open-circuit and 
zero-efficiency conditions. Thus, positive push work is obtained only when the power density exceeds that 
associated with the zero-efficiency condition where j y =Tj a =A n =0 so that 


Portions of the constant-/ 5 circles lying within the positive push work regime are illustrated in 
figure 5. Inspection of equation (40) reveals that the radii of these circles are larger than those that would 
occur for an ideal device in which there is no parasitic voltage drop (A=0). Thus, additional power must be 
applied in a real device to achieve a desired intensity of acceleration. 

It is important to note that the optimum accelerator efficiency occurs with a finite axial current. 
Increasing the accelerator power density until j x vanishes will inevitably reduce the electrical efficiency 
below the optimal value. In fact, to accelerate plasma efficiently (with little heat production), it is desirable 
to maintain gentle acceleration levels by keeping the back emf only slightly less than the applied electric 
field. In practice, the added length and weight associated with gentle acceleration must be traded against 
the electrical inefficiencies encountered with high push power. 

For a generator, the maximum power condition can differ considerably from the maximum effi- 
ciency condition. Furthermore, maximum generator efficiency operation also occurs with a finite axial 
current as observed in figure 4. It should be pointed out that achieving optimal power generation for a 
given value of Hall parameter depends on careful selection of the wall angle. This is examined for the sake 
of completeness. 

By definition, the dimensionless power density may be expressed in the form 


(Q. + (p)j x +{0.(p-l)jy-l + A . 


(47) 



(48) 



( 49 ) 
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If we now use equation (26) to eliminate I in the appropriate expressions for j x ,j y , and E v we obtain the 
dimensionless generator power density in terms of A 


(50) 


) 2 I 

' A “A 0 \ 

o 

<1 

< 

l + A " 

)> 

v Aj — A 0 J 

_ ^ 0 

1 + A 5 _ 


The maximum power density is defined by the constraint dP/d A=0, which yields the relation 


A 0 + A s 


(51) 


Substituting this result into equation (50) then yields 

-(OV) (l + A^l + Aj 


(>v)( 


1 + Q" 


(52) 


Thus, P nvdX depends on (p, and the optimal wall angle for a given Hall parameter is determined by the 
constraint dP max /d(p= 0, which leads to the design criterion, 


Q<p = -1 . 

If this criterion is enforced, we find the maximum possible generator power density to be 


(53) 


P max(£2<p=-l) — 


(l + A lV )(1 + a 0 ) 


(54) 
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6. ELECTRODE BOUNDARY LAYER EFFECTS 


The effective voltage drop arises from cold-wall boundary layer effects including gasdynamic varia- 
tions, discharge constriction, and electrode falls. Because fluid temperature and velocity vary rapidly 
approaching the wall, the conductivity and induced voltage also exhibit strong variations in the thermal 
and momentum boundary layers, respectively. Furthermore, the rapid decrease in temperature near the 
electrode surface leads to breakdown and a sudden switch from diffuse current transport to constricted 
arcs. Current attachment to cold electrodes almost always occurs through a thin layer of short arcs. 

The basic boundary layer features are depicted in figure 6 for an accelerator configuration. For 
simplicity, the thermal and momentum boundary layers are both assumed to have a thickness, 8. In 
addition, current constriction is indicated at the electrode surface where attachment occurs through short 
arcs spanning the anode thickness, 8 a , and the cathode thickness, <5... Note that the illustration shows a 
concentration of current due to the Hall effect at the upstream end of the cathode and the downstream end 
of the anode in accordance with experimental observations for accelerators. 19 ’ 20 

The effective electrode voltage drop, V,, is defined as the potential difference by which the voltage 
between opposite electrodes is increased (due to boundary layer effects) with respect to the fully homoge- 
neous case. The increase in V. is quantified through an integration of the transverse electric field defect 
across the electrode boundary layers. For convenience, V, is split into separate components, 

V d =V a +V c +V g , (55) 

where the anode voltage fall, V ; the cathode voltage fall, V ; and the gasdynamic voltage drop, V , are 

Cl t 5 

defined by 


V a = f 0 ‘(Ey~E y() )dy 

v c = lL c ( E y~ E y«) dy (56) 

V s = . [ E y - E yo) c ( v + j h _ s [ E y “ E y o)dy . 

Note that the zero subscript refers to the homogeneous core flow conditions. 

Factors contributing to the anode and cathode voltage falls include the arc attachment region and 
the potential associated with electron emission from cold electrode material. In general, the voltage drop in 
the arc columns is negligible since the arc length is small compared to the boundary layer thickness, and 
the anode electrode fall is negligible since it does not have to emit electrons. Thus, V a ~0. The cathode does 
emit electrons, however, and the voltage drop associated with this process can only be approximated in an 
empirical fashion. For most materials, V c is approximately 10 to 20 V. 
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The gasdynamic voltage drop may be evaluated by applying the generalized Ohm’s law in the 
boundary layer to eliminate E . In this particular case, the appropriate form for E may be obtained from 
equation (9) by neglecting A: 


^1 + Q 2 ^ 


i\, - HE x + uB . 


Jy 


(57) 


Assuming Q is invariant, V takes the form 


v s=l 


h 


1/1+ n 


+ uB 


0 


r ^ ~8 C 

+ I _ f 

j r, + llB 


J h -8 

X/l + Hq 



7 ~ 2 + “0* 

2o/l + Oo 


J 


>dy 


— ~r~ — y + 

Sq/I + Oo 


>dy 


(58) 


And, it is convenient to split V into two parts 


2 h 



AM 

u a 

[r J 


dy 


(59) 


and 


V a = 2 uqB\ 1 


ff 


j fi 

u c 

«0 


dy 


(60) 


where V 2 and V u account separately for the decrease of X and induced voltage in the boundary layer, and 
X'=X/(1 +Qq 2 ). These integrals may be further simplified by defining %=y/8 and assuming 8 a = 8 C to yield 


y E = ^f' 
1 S'o k 


^-1 

S' 


dt 


(61) 


and 


V u 



(62) 


The lower limit of integration q ( = 8 c /8 is the distance from the electrode surface where the current 
becomes constricted. 
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Evaluation of these integrals is facilitated by the introduction of model boundary layer profiles for 
velocity and temperature. 21 For example, the velocity profile in a turbulent boundary layer is well approxi- 
mated by the power law u/u 0 =(y/8)^^ n , where n~l . If this relation is substituted into the integral for V u and 
£ c =0 is used as the lower limit of integration (since u goes to zero at the wall), we obtain 


V u ~2uoB8/(l + n) . («) 

Evaluation of the integral for is complicated by the fact that X' () /Z' increases steeply near the 
cold electrode surface, which makes the integral extremely sensitive to the lower limit of integration. It is 
assumed that electrical breakdown leading to arc attachment occurs when j /If exceeds the critical electric 
field strength, E c . In this case, the limit, q c , corresponds to a critical conductivity ratio, (0 C , where 


co c 



(64) 


Experience shows that very near the wall the thermal boundary layer may be approximated by the 
power law T/T 0 =(y/8)^ m , where m =18-25. Furthermore, it has also been established that the electrical 
conductivity can be approximated by the simple power law I' /I' q=(T /T^) a {p/p^) b , where a =10-20 and 
b-- 1/2. Thus, for negligible pressure variations, the integral for can be expressed in the simplified form 


■ <65) 

The calculation procedure is as follows. For a given value of j y fl \ 0 and a specified E c , it is possible 
to compute co c , which may be used in turn to compute T ( JT {) =( 0 ( ^ a . Then, the lower limit of integration 
may be determined from the relation £ C =(T C /T 0 ) m and the integral of equation (65) may be evaluated. For 
illustrative purposes, calculated values for q ( . and V^/28 are shown in figure 7 as a function of j\,/I ' 0 , 
assuming E= 40 kV/cm, n=l, m= 18, and <3=10. Note that Vj°^j y in accordance with empirical observations. 
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7. PRACTICAL DESIGN CONSTRAINTS 


In a real device, practical design constraints limit the available operating range. The axial electric 
field, for instance, must be limited due to breakdown considerations. If the axial electric field becomes too 
large, inter-electrode arcing can occur, and the resulting shorts can degrade device performance. Further- 
more, high heat dissipation in these concentrated arc discharges can cause severe erosion of wall material 
and dramatic reduction in channel operating life. 

Experience has shown that a realistic upper limit for the axial electric field is no more than 
100 V/cm. The existence of such a design constraint is illustrated by the dashed lines marked ±E X lim on the 
performance map of figure 5. This implies that DCW accelerators operating with large Q are confined to 
low power-density operation. High power densities are attainable with low Q at the expense of decreased 
efficiency. 

An additional design constraint arises from the observation that excessive heat generation and 
material erosion will occur if the current density entering the electrode surface becomes too large. The 
actual limiting value depends on electrode material and geometry, but experience indicates good channel 
durability for values up to about 10 A/cm 2 . For high power-density applications having short lifetime 
requirements, this value can be exceeded considerably. 
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8. CONCLUSIONS 


The classical cross-plane averaged performance theory for DCW MHD accelerators has been 
extended to include a current-dependent sheath layer voltage drop. This approach yields analytical perfor- 
mance relationships and diagrams that can be used to illuminate the rudimentary behavior of these devices 
and underscore the fundamental interplay of basic parameters. This simplified theoretical treatment is not 
intended for detailed performance prediction of practical devices but rather to aid in developing an intui- 
tive understanding of device operation of general value to analysts, designers, and experimentalists, alike. 
Despite its limitations, the theory can be extremely useful in defining anticipated performance ranges 
while accounting for critical nonideal effects. 
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Table 1. Summary of electrical parameter relationships for diagonalized cross-field MHD devices.* 


Electrical Parameter 

In terms of applied electric field, E x 

In terms of applied electric current, / 

h 

(l-Q(p)l£ x + QZ»5(l + A) 

(l- Q<p)//A + (pZuB( 1+ A) 

1 + Q 2 

l + <p 2 

jy 

(Q + (p)AE x - AuB( 1 + A) 

(Q + <p)//A-2»5(l + A) 

1 + Q 2 

1 + (p 2 

r 

r 

(l + Q 2 )//A + ((p- n)TuB{l+ A) 


x 

2(1 + +) 

E y 

<PE X 

<PE X 

I 

(l + (p 2 )AE x A + (Q- (p)AuB{l + A)A 

I 


1 + Q 2 



*The sign of ip and A depends on the mode of operation. Generator mode: <p <0 and -1<A<0; 
accelerator mode: (p > 0 andA>0. 


Table 2. Summary of mean power relationships for diagonalized cross-field MHD devices.* 


Power Parameter 

In terms of applied electric field, E x 


In terms of applied electric current, / 


P = j E 

ze; 

’ (l + (p 2 ) + (Q- (p) uB{ 1 + A )/E x ’ 

I 2 

( 1 + Q j + (cp — QjZ/tZf( 1 + A) A/I 


1 + Q 2 


AA 2 


1 + (p 2 


A = u ' 0 x B) i 

AE x iiB 

~(Q + (p)-uB(l + A)/E x 


uBI 

~(Q + (p)-AuB(l + A) A//’ 


1 + Q 2 



A 

l + <p 2 


u-(jxB) x 1 

uB 

(Q + (p)-uB{\ + A)/E x 


AuB 


(Q + cp)-AuB(l + A) A/ 1 


j'E % 

E x 

(l + (p 2 ) + (Q-(p)uB{l + A)/E 

X 

I/A 

^1 + Q ^ + ( (p — Q)ZwZi(l + A) A/ 1 


*The sign of (p and A depends on the mode of operation. Generator mode: <p< 0 and-l<A<0; 
accelerator mode: <p>0 andA>0. 


Table 3. Reference conditions for variable A. 


Mode 

A' 

A" 

r 

I" 

Generator 

A 0 

A s 

4=o 

4 

Accelerator 

A_- = A S 

A„=0 

4 

4 
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Figure 1. Composite graphic of alternative cross-field MHD accelerator configurations: (a) Linear Hall acceler 
(b) segmented Faraday accelerator, (c) series connected diagonal accelerator, and (d) DCW accelerator. 
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Figure 2. Orientation of field vectors in a diagonally connected cross-field accelerator. 
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Figure 4. Variation of performance parameters based on representative values 
for flow/field conditions, material functions, and wall angle. Assumes 
that <p=±l, B- 4 T, Z=100, Q=2, u = 3 km/s, A o =-0.1, and A v =-0.02. 
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Figure 5. Performance map for DCW MHD devices with current-dependent voltage drop effect. 
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Figure 6. Illustration of current transport and boundary layer dimensions for a DCW cross-field accelerator. 



32 


Figure 7. Variation of q c and V^/28 as a function of j fZ ' 0 . Assumes £ ( -40 kV/cm, n= 7, m= 18, and a=10. 
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